System and method having optimal, systematic q-ary codes for correcting all asymmetric and symmetric errors of limited magnitude

ABSTRACT

A computer-implemented method and computer program product comprising optimal, systematic q-ary codes for correcting all asymmetric and symmetric errors of limited magnitude are provided.

RELATED APPLICATIONS

This application claims the benefit of priority of U.S. Provisional Application No. 61/358,109, filed on Jun. 24, 2010, the entire contents of which application(s) are incorporated herein by reference.

GOVERNMENT LICENSE RIGHTS

This invention was made with government support under Grant Nos. CCF0701452 and CCF0728810 awarded by National Science Foundation. The government has certain rights in the invention.

FIELD OF THE INVENTION

The present invention relates generally to a systems and methods for correcting errors of limited magnitude, and more particularly, but not exclusively, to computer program products and computer-implemented methods for correcting all asymmetric and symmetric errors of limited magnitude.

BACKGROUND OF THE INVENTION

Classical error control codes have been designed under the assumption of binary symmetric errors, i.e. both 1→0 and 0→1 errors can occur during transmission. Nevertheless, errors in some VLSI and optical systems are asymmetric in nature (5, 8). For example, in VLSI circuits and memories, charges may leak with time but new charges will not be added. Thus, a suitable channel model for such systems is the binary asymmetric channel (Z-channel) where errors are of one dominant type known a priori, say 1→0 errors. In Varshamov et al. (9), Varshamov introduced the q-ary asymmetric channel where the channel's input/output symbols are over the alphabet Q={0, 1, . . . , q−1}. Moreover, such channel has the property that, when a symbol aεQ is transmitted, the corresponding received symbol is in the set {0, 1, . . . , a}, assuming a decreasing error.

Similar to the asymmetric channel is the unidirectional channel; the difference is that the type of error is not known a priori. Asymmetric/unidirectional error control codes have been the subject of many research work (3, 7).

Not until recently has the notion of limited magnitude asymmetric errors been introduced (2); we say that a vector (x_(n-1), x_(n-2), . . . , x₀) over Z_(q) suffers an asymmetric error of maximum magnitude/level l≦q−1 if and only if the corresponding channel output (x′_(n-1), x′_(n-2), . . . , x′₀) is such that x_(i)′ε{x_(i)−l, x_(i)−l+1, . . . , x_(i)}, with x_(i)′εQ. FIG. 1 illustrates the difference between the traditional q-ary asymmetric channel and the q-ary asymmetric channel with level l=1.

In Cassuto et al. (4), an interesting application for this special case of q-ary asymmetric channel was pointed out: multi-level flash memories. Unlike traditional single-level flash memories where each cell stores only one bit, multi-level flash memories achieve higher storage capacities and thus lower manufacturing costs by programming the cells into one of q>2 threshold voltage thereby storing log₂ q bits per cell. Nevertheless, increasing the number of threshold levels imposes an important challenge (6): the voltage difference between states is narrowed since—technically—the voltage window is limited. A natural consequence is that reliability issues such as low data retention and read/write disturbs become more significant (4); errors in such cases are typically in one dominant direction and of limited magnitude.

In Ahlswede et al. (1), the authors introduced codes capable of correcting all asymmetric errors of limited magnitude l (or l-AEC codes for short). However, the proposed codes are non-systematic. A systematic code, where the information symbols are separated from the check symbols, is advantageous over a non-systematic code because, in a systematic code, the data processing and encoding/decoding can be done in parallel. Thus, there is a need in the art for a systematic codes for correcting all asymmetric and symmetric errors of limited magnitude.

SUMMARY OF DISCLOSURE

In one of its aspects, the present disclosure provides systems and methods for correcting all q-ary asymmetric errors of limited magnitude l, where l≦q−2. These systems and methods are advantageous in that the information symbols are separable from the check symbols. Moreover, the systems and methods may use the minimum possible redundancy. An optimal systematic code for correcting all symmetric errors of maximum level l (l-SEC code) is also provided, where

${l \leq \frac{q - 2}{2}},$ based on the similarities between l-SEC and l-AEC codes.

For example, an exemplary method includes a computer-implemented method for systematic encoding of an information vector for use in correction of all asymmetric errors of maximum magnitude l, comprising: providing a q-ary information vector having k information digits, (x_(k-1), x_(k-2), . . . , x₀); computing a, the value of (x_(k-1), x_(k-2), . . . , x₀)mod(l+1) as a vector over (l+1) according to the formula a=y_(k-1)(l+1)^(k-1)+y_(k-2)(l+1)^(k-2)+ . . . +y₀(l+1)⁰, where y_(i)=x_(i) mod(l+1), iε{0, 1, . . . , k−1}; representing a in radix

$\left\lceil \frac{q}{l + 1} \right\rceil$ number with r digits: (a_(r-1), a_(r-2), . . . , a₀), where

${r \geq \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil}};$ computing the check part, c, according to the formula c=(c_(r-1), c_(r-2), . . . , c₀), where c_(i)=(l+1)a_(i), ∀iε{0, 1, . . . , r−1}; and outputting an encoded vector that includes the information digits and check part. In addition, a computer program product is provided, comprising a computer-readable medium having a computer readable program code adapted to be executed to implement the above-method.

Still further, the present disclosure provides a computer-implemented method for systematic decoding of an information vector for use in correction of all asymmetric errors of maximum magnitude l, comprising: providing a q-ary encoded vector that includes an information vector having k received information digits, (x′_(k-1), x′_(k-2), . . . , x′₀) and having received check symbols (c′_(r-1), c′_(r-2), . . . , c′₀), where

${r \geq \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil}};$ recovering the check symbols (c_(r-1), c_(r-2), . . . , c₀) by rounding each received check symbol which is not a multiple of (l+1) upwards to the nearest multiple of (l+1); computing a, the value of

$\left( {\frac{c_{r - 1}}{l + 1},\frac{c_{r - 2}}{l + 1},\ldots\mspace{14mu},\frac{c_{0}}{l + 1}} \right)$ as a vector over

$Z_{\lceil\frac{q}{l + 1}\rceil}$ ${{{as}\mspace{14mu} a} = {{\left( \frac{c_{r - 1}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{r - 1}} + {\left( \frac{c_{r - 2}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{r - 2}} + \ldots + {\left( \frac{c_{0}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{0}}}};$ representing a in radix (l+1) number with k digits as y=(y_(k-1), y_(k-2), . . . , y₀); computing e_(i)=(y_(i)−x′_(i))mod(l+1); and computing the corrected codeword, (x_(k-1), x_(k-2), . . . , x₀), as x_(i)=x′_(i)+e_(i), where i=0, 1, . . . , k−1. The present disclosure also provide analogous correction methods relating to symmetric errors of maximum magnitude l.

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing summary and the following detailed description of the preferred embodiments of the present invention will be best understood when read in conjunction with the appended drawings, in which FIG. 1 illustrates q-ary asymmetric channel (left) vs q-ary asymmetric channel with level 1 (right).

DETAILED DESCRIPTION OF THE DISCLOSURE

The present disclosure provides methods and systems that address the need for error correction of q-ary systems, i.e., those for which q>2. In this regard, in one of its aspects the present invention provides systematic q-ary codes for correcting all asymmetric and symmetric errors of limited magnitude. For example, the present disclosure provides a computer-implemented method and a computer program product for systematic encoding of an information vector for use in correction of all asymmetric errors of maximum magnitude l. The rationale and justification for such methods may be understood in view of the following exposition detailing the formulation of the methods and systems of the present disclosure.

Turning first to a few preliminary principles upon which the present disclosure builds, it turns out that the knowledge of the maximum error level gives nice properties that can be used in the design of error correcting codes. We start by defining a distance metric capturing these properties as introduced in R. Ahlswede et al. (1)

By definition, let x=(x_(n-1), x_(n-2), . . . , x₀) and y=(y_(n-1), y_(n-2), . . . , y₀) be two vectors over Z_(q), then the distance between x and y is defined as: d(x,y)=max{|x _(i) −y _(i) |:iε{0,1, . . . ,n−1}}. It is worth noting that, by properties of absolute values, d(x,y) defines a metric. Moreover, d(x,y)≦q−1 by definition.

The following theorem gives the necessary and sufficient conditions on the minimum distance of an l-AEC code.

Theorem 1 (1): A code C⊂Q^(n) is an l-AEC code if and only if, for all distinct codewords, x, yεC, d(x,y)≧l+1.

Let A_(a)(n,l,q) denote the maximum number of codewords in a q-ary l-AEC code of length n. A bound on A_(a)(n,l,q) and a non-systematic l-AEC code achieving this bound are given in (1):

Theorem 2 (1): ∀nε{1, 2, 3 . . . } and lεQ,

${A_{a}\left( {n,l,q} \right)} = {\left\lceil \frac{q}{l + 1} \right\rceil^{n}.}$ Theorem 3 (1): Let C be a code of length n over Z_(q) defined as C={(x _(n-1) ,x _(n-2) , . . . ,x ₀):x _(i)≡0 mod(l+1), ∀iε{0,1, . . . ,n−1}} Then, C is an l-AEC code with

$\left\lceil \frac{q}{l + 1} \right\rceil^{n}$ codewords.

Finally, for further analysis, x mod a denotes the component-wise remainder of a vector x when divided by an integer a.

Lower Bound on the Number of Check Digits

In the following theorem, the minimum number of check digits needed to encode information vectors of a certain length is presented, which serves, in part, to validate the optimal nature of an l-AEC method/code of the present disclosure.

Theorem 4: Let C be a systematic q-ary l-AEC code, such that the number of information digits in a codeword is k. Then, the number of check digits, r, satisfies the following condition,

$r \geq {\frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil}.}$ Proof. Consider the subset of information vectors, V={(x_(k-1), x_(k-2), . . . , x₀):0≦x_(i)≦l, ∀iε{0, 1, . . . , k−1}}. Vectors of V can be viewed as the set of all vectors of length k over Z_(i+1). Hence, ∀x, yεV, d(x,y)≦l, and |V|=(l+1)^(k). Therefore, by Theorem 1, the checks assigned to vectors in V must be at least l+1 apart for errors to be successfully corrected. Theorem 1 and Theorem 2.2 together give a bound on the number of vectors satisfying such criterion and we get:

$\left\lceil \frac{q}{l + 1} \right\rceil^{r} \geq \left( {l + 1} \right)^{k}$

Taking the log on both sides of the above inequality we get the desired property. Q.E.D.

One important implication of the above theorem is that it is not possible to design a systematic code correcting all errors of maximum magnitude l when l=q−1, since the expression

$\frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil}$ goes to infinity in this case. Therefore, for the following discussion of the l-AEC case, an assumption is made that l≦q−2. An Optimal Systematic l-AEC Code

A method relating to an encoding algorithm is presented for codes which require exactly

$r = \left\lceil \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil} \right\rceil$ check digits, though the present disclosure can also encompasses codes in which

$r \geq {\left\lceil \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil} \right\rceil.}$ l-AEC Encoding Algorithm

Input: The information vector: (x_(k-1), x_(k-2), . . . , x₀)

Output: The encoded vector: x=(x_(k-1), x_(k-2), . . . , x₀, c_(r-1), c_(r-2), . . . , c₀)

A method of the present invention for l-AEC encoding may then include the steps of:

-   -   1. Computing a, the value of (x_(k-1), x_(k-2), . . . ,         x₀)mod(l+1) as a vector over (l+1):         a=y _(k-1)(l+1)^(k-1) +y _(k-2)(l+1)^(k-2) + . . . +y ₀(l+1)⁰,         where         -   y_(i)=x_(i) mod(l+1), iε{0, 1, . . . , k−1}.         -   That is, a is the number with radix (l+1) representation of             (x_(k-1), x_(k-2), . . . , x₀)mod(l+1).     -   2. Representing a in radix

$\left\lceil \frac{q}{l + 1} \right\rceil$ number with r digits: (a_(r-1), a_(r-2), . . . , a₀).

-   -   3. Computing the check part:         c=(c _(r-1) ,c _(r-2) , . . . ,c ₀), where c _(i)=(l+1)a _(i) ,         ∀iε{0,1, . . . ,r−1}.     -   4. Outputting the encoded vector x=(x_(k-1), x_(k-2), . . . ,         x₀, c_(r-1), c_(r-2), . . . , c₀).

EXAMPLE 1

The word (6,2,8,1) is encoded over Z₁₀, assuming a maximum error level of 2. The minimum number of check digits needed is

$r = {\left\lceil \frac{4{\log\left( {2 + 1} \right)}}{\log\left\lceil \frac{10}{2 + 1} \right\rceil} \right\rceil = 4.}$ With notations as above, a=0×3³+2×3²+2×3¹+1×3⁰=25 and thus (0,1,2,1) is the representation of a in base 4. Therefore, the encoded codeword is (6,2,8,1,0,3,6,3).

The following theorem further supports the optimal nature of the above method. Theorem 5: The above construction yields codewords of minimum distance l+1. Proof. Given two distinct information vectors: ν₁=(x_(k-1), x_(k-2), . . . , x₀) and ν₂=(y_(k-1), y_(k-2), . . . , y₀), there are two possibilities: ν₁ mod(l+1)≠ν₂ mod(l+1) in which case each vector is assigned different check digits and, since the check digits are multiples of (l+1), the distance between the resulting codewords is at least l+1. For the second case, ν₁ mod(l+1)=ν₂ mod(l+1), ν₁ and ν₂ are assigned the same check digits. Nevertheless, by distinctness of ν₁ and ν₂, ∃iε{0, 1, . . . , k−1} such that x_(i)=y_(i)+m(l+1), m≧1. Therefore, d(ν₁,ν₂)≧l+1, and the resulting codewords satisfy the desired property. Q.E.D.

l-AEC Decoding Algorithm

Input: The channel output: x′=(x′_(k-1), x′_(k-2), . . . , x′₀, c′_(r-1), c′_(r-2), . . . , c′₀)

Output: The recovered codeword: (x_(k-1), x_(k-2), . . . , x₀, c_(r-1), c_(r-2), . . . , c₀)

A method of the present invention for l-AEC decoding may then include the steps of:

-   -   1. Recovering the check symbols, (c_(r-1), c_(r-2), . . . , c₀),         by rounding each received check symbol which is not a multiple         of (l+1) upwards to the nearest multiple of (l+1).     -   2. Computing a, the value of

$\left( {\frac{c_{r - 1}}{l + 1},\frac{c_{r - 2}}{l + 1},\ldots\mspace{14mu},\frac{c_{0}}{l + 1}} \right)$ as a vector over

${Z_{\lceil\frac{q}{l + 1}\rceil}:a} = {{\left( \frac{c_{r - 1}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{r - 1}} + {\left( \frac{c_{r - 2}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{r - 2}} + \ldots + {\left( \frac{c_{0}}{l + 1} \right){\left\lceil \frac{q}{l + 1} \right\rceil^{0}.}}}$

-   -   -   That is, a is the number with radix

$\left\lceil \frac{q}{l + 1} \right\rceil$ representation of

$\left( {\frac{c_{r - 1}}{l + 1},\frac{c_{r - 2}}{l + 1},{\ldots\mspace{14mu}\frac{c_{0}}{l + 1}}} \right).$

-   -   3. Representing a in radix (l+1) number with k digits:         y=(y _(k-1) ,y _(k-2) , . . . , y ₀).     -   4. Letting e_(i)=(y_(i)−x′_(i))mod(l+1), where the information         symbols are::         (x _(k-1) ,x _(k-2) , . . . ,x ₀), such that x _(i) =x′ _(i) +e         _(i) , i=0,1, . . . ,k−1.     -   5. Outputting the codeword: (x_(k-1), x_(k-2), . . . , x₀,         c_(r-1), c_(r-2), . . . , c₀).

EXAMPLE 2

Let the encoded word be as in Example 1 and the channel output be x′=(4,2,7,1,0,3,5,1). Rounding the check symbols which are not multiple of three upwards to the nearest multiple of 3, we get (0,3,6,3). As in Steps 2 and 3 of the algorithm, we compute a=0121₄=25, and y=(0,2,2,1). Thus, the correct information symbols are (6,2,8,1).

The following theorem validates the above decoding method. Theorem 6: Let x be a codeword encoded using the encoding algorithm given in above, and let x′ be the l-asymmetric channel output. Then, the above l-AEC decoding method successfully recovers x. Proof Let x=(x_(k-1), x_(k-2), . . . , x₀, c_(r-1), c_(r-2), . . . , c₀), then, by the channel properties, x′=(x′_(k-1), x′_(k-2), . . . , x′₀, c′_(r-1), c′_(r-2), . . . , c′₀) is such that x_(i)−l≦x′_(i)≦x_(i) and c_(i)−l≦c′_(i)≦c_(i). Moreover, the encoding algorithm yields check symbols that are multiples of l+1, i.e. c′_(i)lies between two successive multiples of l+1. Therefore, the first step of the decoding algorithm successfully recovers the check symbols. It can also be seen that Steps 2 and 3 of the above decoding algorithm do the reverse operations of Steps 1 and 2 of the encoding algorithm. Hence, with notations as above, y=(y_(k-1), y_(k-2), . . . , y₀)=(x_(k-1), x_(k-2), . . . , x₀)mod(l+1). At Step 4, e_(i) can be seen as the magnitude of the error at the information symbol, such that x _(i) =x′ _(i) +e _(i) ≡y _(i)(mod(l+1)).

Thus, e≡(y _(i) −x′ _(i))(mod(l+1)). Since the maximum error magnitude is l(i.e. 0≦e_(i)≦l) the value of e_(i) is successfully computed at Step 4 of the decoding algorithm, recovering the information symbols. Q.E.D. An Optimal Systematic l-SEC Code

In addition to methods relating to correction of all asymmetric errors of limited magnitude, the present invention provides methods for correcting symmetric error of l limited magnitude.

For a vector x=(x_(n-1), x_(n-2), . . . , x₀) over Z, the corresponding channel output x′=(x′_(n-1), x′_(n-2), . . . , x′₀) suffers a symmetric error of maximum magnitude l if and only if x_(i)−e_(i)≦x′_(i)≦x_(i)+e_(i), where 0≦e_(i)≦l, ∀iε{0, 1, . . . , n}. As shown below, the similarity between the properties of l-SEC codes and l-AEC codes allows extension of the l-AEC methods to a family of l-SEC methods.

The following theorems serve to set conditions that should be met for an optimal l-SEC code. Theorem 7: A code C is capable of correcting all symmetric errors of maximum magnitude l if and only if C has minimum distance 2l+1. Proof. Let S_(x) be the set of all words obtained from a codeword xεC, where |x|=n, due to n or less symmetric errors of maximum magnitude l, i.e. S_(x)={(x′_(n-1), x′_(n-2), . . . , x′₀):x′_(i)=x_(i)±e_(i), e_(i)ε{0, 1 . . . , l}}. Then, it is easy to see that, if C has minimum distance 2l+1, then ∀x, yεC, S_(x)∩S_(y)=Ø. Therefore, C is an l-SEC code. Conversely, if ∃x, yεC such that d(x,y)≦2l, then it is possible to obtain a word, say z, from both x and y due to symmetric errors of magnitude l or less. Hence, a decoder for C cannot correct z. Therefore, the minimum distance should be no less than 2l+1 for C to be capable of correcting all symmetric errors of maximum magnitude l. Q.E.D.

The above theorem implies that, when

${l > \frac{q - 2}{2}},$ any l-SEC code can have at most one codeword. Thus, we assume

$l \leq {\frac{q - 2}{2}.}$

Theorem 8: Let A (n,l,q) denote the maximum number of words in a q-ary l-SEC code of length n. Then

${A_{s}\left( {n,l,q} \right)} = {\left\lceil \frac{q}{{2l} + 1} \right\rceil^{n}.}$ Proof Similar to the proof of Theorem 2. Q.E.D.

Theorem 9: Let C be a code of length n over Z_(q) defined as: C={(x _(n-1) ,x _(n-2) , . . . ,x ₀):x _(i)≡0(mod(2l+1), ∀iε{0,1, . . . ,n−1}}.

Then, C is an l-SEC with A_(s)(n,l,q) codewords. Proof. Similar to the proof of Theorem 3. Digits of the channel output can be decoded in this case by rounding downwards or upwards (whichever is closer) to the nearest multiple of 2l+1. Q.E.D.

Theorem 10: Let C be a systematic q-ary l-SEC code, such that the number of information digits in a codeword is k. Then, the number of check digits, r, satisfies the following condition

$r \geq {\frac{k \times {\log\left( {{2l} + 1} \right)}}{\log\left\lceil \frac{q}{{2l} + 1} \right\rceil}.}$

Proof. Similar to the proof of Theorem 4. Q.E.D.

Now that we have identified the similarities between l-SEC and l-AEC codes, it can easily be seen that the proposed encoding l decoding algorithms for l-AEC codes can be modified in the following way to construct an l-SEC code.

l-SEC Encoding Algorithm

Input: The information vector: (x_(k-1), x_(k-2), . . . , x₀)

Output: The encoded vector: x=(x_(k-1), x_(k-2), . . . , x₀, c_(r-1), c_(r-2), . . . , c₀)

A method of the present invention for l-SEC encoding may then include the steps of:

-   -   1. Computing a, the value of (x_(k-1), x_(k-2), . . . ,         x₀)mod(2l+1) as a vector over (2l+1):         a=y _(k-1)(2l+1)^(k-1) +y _(k-2)(2l+1)^(k-2) + . . . +y         ₀(2l+1)⁰,         -   where y_(i)=x_(i) mod(2l+1), iε{0, 1, . . . , k−1}     -   2. Representing a in radix

$\left\lceil \frac{q}{{2l} + 1} \right\rceil$ number with r digits: (a _(r-1) ,a _(r-2) , . . . ,a ₀).

-   -   3. Computing the check part:         c=(c _(r-1) ,c _(r-2) , . . . ,c ₀).         -   where c_(i)=(2l+1)a_(i), ∀iε{0, 1, . . . , r−1}.     -   4. Outputting the encoded vector x=(x_(k-1), x_(k-2), . . . ,         x₀, c_(r-1), c_(r-2), . . . , c₀)

EXAMPLE 3

We encode the word (5,3,7,0) over Z₁₀, assuming a maximum symmetric error level of 1. The number of check digits needed is

$r = {\left\lceil \frac{4{\log\left( {{2 \times 1} + 1} \right)}}{\log\left\lceil \frac{10}{{2 \times 1} + 1} \right\rceil} \right\rceil = 4.}$ Then, (y₃, y₂, y₁, y₀)=(2,0,1,0) and a=2×3³+0×3²+1×3¹+0×3⁰=57. Now, a=57 in radix

$\left\lceil \frac{q}{{2l} + 1} \right\rceil = {\left\lceil \frac{10}{4} \right\rceil = 4}$ representation is (0,3,2,1) and so, the check is (0,9,6,3). Therefore, the encoded codeword is (5,3,7,0,0,9,6,3). l-SEC Decoding Algorithm

Input: The channel output: x′=(x′_(k-1), x′_(k-2), . . . , x′₀, c′_(r-1), c′_(r-2), . . . , c′₀)

Output: The recovered codeword: (x_(k-1), x_(k-2), . . . , x₀, c_(r-1), c_(r-2), . . . , c₀)

A method of the present invention for l-SEC decoding may then include the steps of:

-   -   1. Recovering the check symbols, (c_(r-1), c_(r-2), . . . , c₀),         by rounding each received check symbol which is not a multiple         of 2l+1 either upwards or downwards to the nearest multiple of         2l+1, whichever is closer.     -   2. Computing a, the value of

$\left( {\frac{c_{r - 1}}{{2l} + 1},\frac{c_{r - 2}}{{2l} + 1},\ldots\mspace{14mu},\frac{c_{0}}{{2l} + 1}} \right)$ as a vector over

${Z_{\lceil\frac{q}{{2l} + 1}\rceil}:a} = {{\left( \frac{c_{r - 1}}{{2l} + 1} \right)\left\lceil \frac{q}{{2l} + 1} \right\rceil^{r - 1}} + {\left( \frac{c_{r - 2}}{{2l} + 1} \right)\left\lceil \frac{q}{{2l} + 1} \right\rceil^{r - 2}} + \ldots + {\left( \frac{c_{0}}{{2l} + 1} \right)\left\lceil \frac{q}{{2l} + 1} \right\rceil^{0}}}$

-   -   3. Representing a in radix (2l+1) number with k digits:         y=(y _(k-1) ,y _(k-2) , . . . ,y ₀), where the information         symbols are         (x _(k-1) ,x _(k-2) , . . . ,x ₀), such that ∀iε{0,1, . . .         ,k−1} x _(i) =x′ _(i) +e _(i), where         e _(i)≡(y _(i) −x′ _(i))mod(2l+1) such that −l≦e _(i) ≦l and         x _(i) ≡y _(i)(mod(2l+1)).     -   5. Outputting the codeword: (x_(k-1), x_(k-2), . . . , x₀,         c_(r-1), c_(r-2), . . . , c₀).

EXAMPLE 4

Let the encoded word be as in Example 3 and the channel output be x′=(4,2,8,1,1,8,7,2). Rounding the check symbols (that are not multiple of 3) upwards or downwards to the nearest multiple of 3 we get (0,9,6,3). As in Steps 2 and 3 of the l-SEC algorithm, we compute a=0321₄=57, and 57 in radix (2l+1)=3 representation is y=(y₃,y₂,y₁,y₀)=(2,0,1,0). Thus the error vector is e=(e ₃ ,e ₂ ,e ₁ ,e ₀)=(y ₃ −x′ ₃ ,y ₂ −x′ ₂ ,y ₁ −x′ ₁ ,y ₀ −x′ ₀)(mod(2l+1))=((2−4)mod3,(0−2)mod3,(1−8)mod3,(0−1)mod3)=(1,1,−1,−1). Then the correct information word is (4+1,2+1,8−1,1−1)=(5,3,7,0).

Since our code constructions and resulting methods are systematic, information symbols are separable from the check symbols resulting in faster encoding l decoding operations. Fortunately, the cost of having systematic code is low: we show that the rate of our l-AEC code is very close to the one given in (1). The rate R of an error correcting code is given as

$R = \frac{k}{n}$ where k is the number of information digits and n is the length of the code. For the code given in (1), the rate, R_(non-systematic), is

$\begin{matrix} {R_{{non}\text{-}{systematic}} = \frac{\log_{q}\left\lceil \frac{q}{l + 1} \right\rceil^{n}}{n}} \\ {= {\log_{q}{\left\lceil \frac{q}{l + 1} \right\rceil.}}} \end{matrix}$

The rate of the proposed code, R_(systematic), is

$\begin{matrix} {R_{systematic} = \frac{k}{k + r}} \\ {= \frac{k}{k + {k\left\lceil \frac{\log_{q}\left( {l + 1} \right)}{\log_{q}\left\lceil \frac{q}{l + 1} \right\rceil} \right\rceil}}} \\ {{= \frac{1}{1 + \left\lceil \frac{\log_{q}\left( {l + 1} \right)}{\log_{q}\left\lceil \frac{q}{l + 1} \right\rceil} \right\rceil}},} \end{matrix}$

which can be approximated to

$\begin{matrix} {R_{systematic} \approx \frac{1}{1 + \frac{\log_{q}\left( {l + 1} \right)}{\log_{q}\frac{q}{l + 1}}}} \\ {= \frac{\log_{q}\frac{q}{l + 1}}{{\log_{q}\frac{q}{l + 1}} + {\log_{q}\left( {l + 1} \right)}}} \\ {= \frac{\log_{q}\frac{q}{l + 1}}{{\log_{q}q} - {\log_{q}\left( {l + 1} \right)} + {\log_{q}\left( {l + 1} \right)}}} \\ {= {\log_{q}{\frac{q}{l + 1}.}}} \end{matrix}$

The table below illustrates the closeness of the values of R_(systematic) and R_(non-systematic) for q=8 and q=16:

q l R_(non-systematic) R_(systematic) 8 1 0.66 0.5 2 0.52 0.5 3 0.33 0.33 4 0.33 0.25 5 0.33 0.25 6 0.33 0.25 16 1 0.75 0.5 2 0.64 0.5 3 0.5 0.5 4 0.5 0.33 5 0.39 0.33 6 0.39 0.33 7 0.25 0.25 8 0.25 0.25 9 0.25 0.25 10 0.25 0.25 11 0.25 0.25 12 0.25 0.25 13 0.25 0.25 14 0.25 0.25

These and other advantages of the present invention will be apparent to those skilled in the art from the foregoing specification. Accordingly, it will be recognized by those skilled in the art that changes or modifications may be made to the above-described embodiments without departing from the broad inventive concepts of the invention. It should therefore be understood that this invention is not limited to the particular embodiments described herein, but is intended to include all changes and modifications that are within the scope and spirit of the invention as set forth in the claims.

Several publications are referenced in this application, the disclosures of each of these publications are incorporated by reference herein.

REFERENCES

(1) R. Ahlswede, H. Aydinian, L. H. Khachatrian, and L. M. G. M. Tolhuizen, “On q-ary codes correcting all unidirectional errors of a limited magnitude”, Proceedings of Ninth International Workshop on Algebraic and Combinatorial Coding Theory, Kranevo, Bulgaria, pp. 20-26, June 2004.

(2) R. Ahlswede, H. Aydinian and L. H. Khachatrian, “Unidirectional error control codes and related combinatorial problems”, Proceedings of Eight International Workshop on Algebraic and Combinatorial Coding Theory, Russia, pp. 6-9, September 2002.

(3) M. Blaum, “Codes for Detecting and Correcting Unidirectional Errors”, IEEE Computer Society Press, Los Alamitos, Calif., 1993.

(4) Y. Cassuto, M. Schwartz, V. Bohossian and J. Bruck, “Codes for multi-level flash memories: correcting asymmetric limited-magnitude errors”, Proceedings of ISIT2007, Nice, France, June 2007.

(5) S. D. Constantin and T. R. N. Rao, “On the theory of binary asymmetric error correcting codes”, Information and computation, vol. 40, pp. 20-36, 1979.

(6) B. Eitan, R. Kazerounian, A. Roy, G. Crisenza, P. Cappelletti and A. Modelli, “Multilevel flash cells and their trade-offs”, International Electron Devices Meeting, pp. 169-172, December 1996.

(7) T. Klove, “Error correcting codes for the asymmetric channel”, In Dept. Mathematics, Univ. Bergen, Bergen, Norway, Tech. Rep. 1809-0781, 1995.

(8) J. R. Pierce, “Optical Channels: Practical Limits with Photon Counting”, IEEE Trans. Communications, vol. 26, no. 12, pp. 1819-1821, December 1978.

(9) R. R. Varshamov, “A class of codes for asymmetric channels and a problem from the additive theory of numbers”, IEEE Trans. Inform. Theory, vol. 19, no. 1, pp. 92-95, January 1973. 

What is claimed is:
 1. A computer-implemented method for systematic encoding of an information vector for use in correction of all asymmetric errors of maximum magnitude l and executed by a processor, comprising: a) providing a q-ary information vector having k information digits, (x_(k-1),x_(k-2), . . . ,x₀); b) computing a, the value of (x_(k-1),x_(k-2), . . . ,x₀) mod (l+1) as a vector over (l+1) according to the formula a=y _(k-1)(l+1)^(k-1) +y _(k-2)(l+1)^(k-2) +. . . +y ₀(l+1)⁰, where y _(i)=x _(i)mod(l+1),iε{0,1, . . . ,k−1}; c) representing a in radix $\left\lceil \frac{q}{l + 1} \right\rceil$  number with r digits: (a_(r-1),a_(r-2), . . . ,a₀), where ${r \geq \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil}};$ d) computing the check part, c, according to the formula c=(c_(r-1),c_(r-2), . . . ,c₀), where c_(i)=(l+1)a_(i),∀iε{0,1, . . . , r−1}; and e) outputting an encoded vector that includes the information digits and check part, wherein any two check parts computed according to step (d) are at least l+1 apart.
 2. The method according to claim 1, wherein $r = {\left\lceil \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil} \right\rceil.}$
 3. A computer program product, comprising a non-transitory computer-readable medium having a computer readable program code embodied therein, said computer readable program code adapted to be executed to implement a method for encoding an information vector for use in systematic correction of all asymmetric errors of maximum magnitude /comprising:A=πr² a) providing a q-ary information vector having k information digits, (x_(k-1),x_(k-2), . . . ,x₀); b) computing a , the value of (x_(k-1),x_(k-2), . . . ,x₀) mod (l+1) as a vector over (l+1) according to the formula a=y _(k-1)(l+1)^(k-1) +y _(k-2)(l+1)^(k-2) +. . . +y ₀(l+1)⁰, where y _(i)=x _(i)mod(l+1),iε{0,1, . . . ,k−1}; c) representing a in radix $\left\lceil \frac{q}{l + 1} \right\rceil$  number with r digits: (a_(r-1),a_(r-2), . . . ,a₀), where ${r \geq \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil}};$ d) computing the check part, c, according to the formula c=(c_(r-1),c_(r-2), . . . ,c₀), where c_(i)=(l+1)a_(i),∀iε{0,1, . . . ,r−1}; and e) outputting an encoded vector that includes the information digits and check part, wherein any two check parts computed according to step (d) are at least l+1 apart.
 4. The computer program product according to claim 3, wherein $r = {\left\lceil \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil} \right\rceil.}$
 5. A computer-implemented method for systematic decoding of an information vector for use in correction of all asymmetric errors of maximum magnitude / and executed by a processor, comprising: a) providing a q-ary encoded vector that includes an information vector having k received information digits, (x′_(k-1),x′_(k-2), . . . ,x′₀) and having received check symbols (c′_(r-1),c′_(r-2), . . . ,c′₀), where ${r \geq \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil}};$ b) recovering the check symbols (c_(r-1),c_(r-2), . . . ,c₀) by rounding each received check symbol which is not a multiple of (l+1) upwards to the nearest multiple of (l+1); c) computing a , the value of $\left( {\frac{c_{r - 1}}{l + 1},\frac{c_{r - 2}}{l + 1},\ldots\mspace{14mu},\frac{c_{0}}{l + 1}} \right)$  as a vector over $Z_{\lceil\frac{q}{l + 1}\rceil}$  as ${a = {{\left( \frac{c_{r - 1}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{r - 1}} + {\left( \frac{c_{r - 2}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{r - 2}} + \ldots + {\left( \frac{c_{0}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{0}}}};$ d) representing a in radix (l+1) number with k digits as y=(y_(k-1),y_(k-2) , . . . ,y₀); e) computing e_(i)=(y_(i)−x′_(i)) mod (l+1) ; and computing the corrected codeword, (x_(k-1),x_(k-2), . . . ,x₀), as x_(i)=x′_(i)+e_(i), where i=0,1, . . . ,k−1, wherein any two of the check symbols are at least l+1 apart.
 6. The method according to claim 5, wherein $r = {\left\lceil \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil} \right\rceil.}$
 7. A computer program product, comprising a non-transitory computer-readable medium having a computer readable program code embodied therein, said computer readable program code adapted to be executed to implement a method for decoding an information vector for use in systematic correction of all asymmetric errors of maximum magnitude /comprising: a) providing a q-ary encoded vector that includes an information vector having k received information digits, (x′_(k-1),x′_(k-2), . . . ,x′₀) and having received check symbols (c′_(r-1),c′_(r-2), . . . ,c′₀), where ${r \geq \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil}};$ b) recovering the check symbols (c_(r-1),c_(r-2), . . . ,c₀) by rounding each received check symbol which is not a multiple of (l+1) upwards to the nearest multiple of (l+1); c) computing a , the value of $\left( {\frac{c_{r - 1}}{l + 1},\frac{c_{r - 2}}{l + 1},\ldots\mspace{14mu},\frac{c_{0}}{l + 1}} \right)$  as a vector over $Z_{\lceil\frac{q}{l + 1}\rceil}$  as ${a = {{\left( \frac{c_{r - 1}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{r - 1}} + {\left( \frac{c_{r - 2}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{r - 2}} + \ldots + {\left( \frac{c_{0}}{l + 1} \right)\left\lceil \frac{q}{l + 1} \right\rceil^{0}}}};$ d) representing a in radix (l+1) number with k digits as y=(y_(k-1),y_(k-2), . . . , y₀); e) computing e_(i)=(y_(i)−x′_(i)) mod (l+1) ; and f) computing the corrected codeword, (x_(k-1),x_(k-2), . . . ,x₀), as x′_(i)+e_(i), where i=0,1, . . . ,k−1, wherein any two of the check symbols are at least l+1 apart.
 8. The computer program product according to claim 7, wherein $r = {\left\lceil \frac{k \times {\log\left( {l + 1} \right)}}{\log\left\lceil \frac{q}{l + 1} \right\rceil} \right\rceil.}$
 9. A computer-implemented method for systematic encoding of an information vector for use in correction of all symmetric errors of maximum magnitude / and executed by a processor, comprising: a) providing a q-ary information vector having k information digits, (x_(k-1),x_(k-2), . . . ,x₀); b) computing a , the value of (x_(k-1),x_(k-2), . . . ,x₀)mod(2l+1) as a vector over (2l+1) according to the formula a=y _(k-1)(2l+1)^(k-1) +y _(k-2)(2l+1)^(k-2)+. . .+y ₀(2l+1)⁰, where y_(i)=x_(i)mod(2l+1),iε{0,1, . . . ,k−1}; c) representing a in radix $\left\lceil \frac{q}{{2l} + 1} \right\rceil$  number with r digits, (a_(r-1),a_(r-2), . . . ,a₀), where ${r \geq \frac{k \times {\log\left( {{2l} + 1} \right)}}{\log\left\lceil \frac{q}{{2l} + 1} \right\rceil}};$ d) computing the check part as c =(c_(r-1),c_(r-2), . . . ,c₀), where c_(i)=(2l+1)a_(i),∀iε{0,1, . . . ,r−1}; and e) outputting an encoded vector that includes the information digits and check part, wherein any two check parts computed according to step (d) are at least l+1 apart.
 10. The method according to claim 9, wherein $r = {\frac{k \times {\log\left( {{2l} + 1} \right)}}{\log\left\lceil \frac{q}{{2l} + 1} \right\rceil}.}$
 11. A computer program product, comprising a non-transitory computer-readable medium having a computer readable program code embodied therein, said computer readable program code adapted to be executed to implement a method for encoding an information vector for use in systematic correction of all symmetric errors of maximum magnitude /comprising: a) providing a q-ary information vector having k information digits, (x_(k-1),x_(k-2), . . . ,x₀); b) computing a , the value of (x_(k-1),x_(k-2), . . . ,x₀)mod(2l+1) as a vector over (2l+1) according to the formula a=y _(k-1)(2l+1)^(k-1) +y _(k-2)(2l+1)^(k-2)+. . .+y ₀(2l+1)⁰, where y_(i)=x_(i)mod(2l+1),iε{0,1, . . . ,k−1}; c) representing a in radix $\left\lceil \frac{q}{{2l} + 1} \right\rceil$  number with r digits, (a_(r-1),a_(r-2), . . . ,a₀), where ${r \geq \frac{k \times {\log\left( {{2l} + 1} \right)}}{\log\left\lceil \frac{q}{{2\; l} + 1} \right\rceil}};$ d) computing the check part as c=(c_(r-1),c_(r-2), . . . ,c₀), where c_(i)=(2l+1)a_(i),∀iε{0,1, . . . ,r−1}; and e) outputting an encoded vector that includes the information digits and check part, wherein any two check parts computed according to step (d) are at least l+1 apart.
 12. The computer program product according to claim 11, wherein $r = {\frac{k \times {\log\left( {{2l} + 1} \right)}}{\log\left\lceil \frac{q}{{2l} + 1} \right\rceil}.}$
 13. A computer-implemented method for systematic decoding of an information vector for use in correction of all symmetric errors of maximum magnitude l and executed by a processor, comprising: a) providing a q-ary encoded vector that includes an information vector having k received information digits, (x′_(k-1),x′_(k-2), . . . ,x′₀) and having received check symbols (c′_(r-1),c′_(r-2), . . . ,c′₀), where ${r \geq \frac{k \times {\log\left( {{2l} + 1} \right)}}{\log\left\lceil \frac{q}{{2l} + 1} \right\rceil}};$ b) recovering the check symbols by rounding each received check symbol (c_(r-1),c_(r-2), . . . ,c₀)which is not a multiple of (2l+1) either upwards or downwards to the nearest multiple of (2l+1) , whichever is closer; c) computing a , the value of $\left( {\frac{c_{r - 1}}{{2l} + 1},\frac{c_{r - 2}}{{2l} + 1},\ldots\mspace{14mu},\frac{c_{0}}{{2l} + 1}} \right)$  as a vector over $Z_{\lceil\frac{q}{{2l} + 1}\rceil}$  as ${a = {{\left( \frac{c_{r - 1}}{{2l} + 1} \right)\left\lceil \frac{q}{{2l} + 1} \right\rceil^{r - 1}} + {\left( \frac{c_{r - 2}}{{2l} + 1} \right)\left\lceil \frac{q}{{2l} + 1} \right\rceil^{r - 2}} + \ldots + {\left( \frac{c_{0}}{{2l} + 1} \right)\left\lceil \frac{q}{{2l} + 1} \right\rceil^{0}}}};$ d) representing a in radix (2l+1) number with k digits as y=(y_(k-1),y_(k-2), . . . ,y₀); e) computing e_(i)=(y_(i)−x′_(i))mod(2l+1) such that −l≦e_(i)≦l; and f) computing the corrected codeword, (x_(k-1),x_(k-2), . . . ,x₀), as x_(i)=x′_(i)+e_(i), where i=0,1, . . . ,k−1, wherein any two check symbols are at least l+1 apart.
 14. The method according to claim 13, wherein $r = {\frac{k \times {\log\left( {{2l} + 1} \right)}}{\log\left\lceil \frac{q}{{2l} + 1} \right\rceil}.}$
 15. A computer program product, comprising a non-transitory computer-readable medium having a computer readable program code embodied therein, said computer readable program code adapted to be executed to implement a method for decoding an information vector for use in systematic correction of all symmetric errors of maximum magnitude /comprising: a) providing a q-ary encoded vector that includes an information vector having k received information digits, (x′_(k-1),x′_(k-2), . . . ,x′₀) and having received check symbols (c′_(r-1),c′_(r-2), . . . ,c′₀), where ${r \geq \frac{k \times {\log\left( {{2l} + 1} \right)}}{\log\left\lceil \frac{q}{{2l} + 1} \right\rceil}};$ b) recovering the check symbols (c_(r-1),c_(r-2), . . . ,c₀) by rounding each received check symbol which is not a multiple of (2l+1) either upwards or downwards to the nearest multiple of (2l+1) , whichever is closer; c) computing a , the value of $\left( {\frac{c_{r - 1}}{{2l} + 1},\frac{c_{r - 2}}{{2l} + 1},\ldots\mspace{14mu},\frac{c_{0}}{{2l} + 1}} \right)$  as a vector over $Z_{\lceil\frac{q}{{2l} + 1}\rceil}$  as ${a = {{\left( \frac{c_{r - 1}}{{2l} + 1} \right)\left\lceil \frac{q}{{2l} + 1} \right\rceil^{r - 1}} + {\left( \frac{c_{r - 2}}{{2l} + 1} \right)\left\lceil \frac{q}{{2l} + 1} \right\rceil^{r - 2}} + \ldots + {\left( \frac{c_{0}}{{2l} + 1} \right)\left\lceil \frac{q}{{2l} + 1} \right\rceil^{0}}}};$ d) representing a in radix (2l+1) number with k digits as y=(y_(k-1),y_(k-2), . . . ,y₀); e) computing e_(i)=(y_(i)−x′_(i))mod(2l+1) such that −l≦e_(i)≦l; and f) computing the corrected codeword, (x_(k-1),x_(k-2), . . . ,x₀), as x_(i)=x′_(i)+e_(i), where i=0,1, . . . ,k−1, wherein any two check symbols are at least l+1 apart.
 16. The computer program product according to claim 15, wherein $r = {\frac{k \times {\log\left( {{2l} + 1} \right)}}{\log\left\lceil \frac{q}{{2l} + 1} \right\rceil}.}$ 